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Abstract 

The antiferromagnetic critical point of the Potts model on the square lattice was identified by Baxter [l] as a 
staggered integrable six- vertex model. In this work, we investigate the integrable structure of this model. It enables 
us to derive some new properties, such as the Hamiltonian limit of the model, an equivalent vertex model, and the 
structure resulting from the Zi symmetry. Using this material, we discuss the low-energy spectrum, and relate it to 
geometrical excitations. We also compute the critical exponents by solving the Bethe equations for a large lattice 
width N. The results confirm that the low-energy spectrum is a collection of continua with typical exponent gaps 
of order (log N)~ 2 . 

1 Introduction 

Our understanding of 1 + 1 conformal field theories with non compact target spaces has improved a great deal in the 
last few years, thanks to the use of geometrical methods [2], and ideas from string theory [3J. The topic is of the 
highest interest in the context of the AdS/CFT duality. 

Theories with non compact target spaces also play an important role in statistical mechanics. A sophisticated 
example of this role occurs in the supersymmetric approach to phase transitions for non interacting disordered electronic 
systems, where the universality class of the transition between plateaux of the integer quantum Hall effect is related 
with the IR limit of a non compact 1 + 1 supersigma model at 6 = it [I]. A more basic example is provided by 
Brownian motion and subtle properties thereof, such as the (non) intersection exponents [5]. In both cases, the non 
compacity of the target space occurs because the electron trajectories or the random path can visit a given site (edge) 
an infinite number of times. This is in sharp contrast with self avoiding models for which almost everything is by now 
understood, and related with ordinary CFTs (essentially, a twisted free boson). 

An obvious strategy to tackle the physics of models with non compact target spaces is to start with a lattice model 
having an infinity of degrees of freedom per site/link. For instance, it is easy to generalize the usual XXX chain to a 
non compact representation of SL(2,R), and try to use the standard tools of Bethe ansatz, Baxter Q-operator, etc, 
to obtain properties such as gaplessness and critical exponents. Despite some serious progress in this direction [6] , the 
problem is far from being closed. 
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Another strategy is based on the observation that a non compact continuum limit may well arise from a lattice 
model with finite number of degrees of freedom per site/link if the non unitarity is strong enough. The two families of 
examples exhibited so far involve models with supergroup symmetries — either models with OSP(m\2n) symmetries 
(such as intersecting loop models) [7] which, in their Goldstone phases can be described in the IR by a collection of 
free bosons and symplectic fermions, or the SU(l\2) integrable spin chain with alternating 3 and 3 representations, 
which was found to be described by the SU (2\1) WZW model at k = 1 [5]. In both cases, a continuous spectrum of 
critical exponents is found, and the target space does exhibit some non compact directions indeed. 

The examples of Brownian motion and self intersecting dense curves should convince the reader that non compact 
target spaces are more common and useful than might have been surmised a few years ago. In a recent paper, we found 
[9] that the antiferromagnetic Potts model on the square lattice for Q continuous has critical properties seemingly 
involving a twisted non compact boson. This conclusion was based on some numerical and analytical evidence, and 
implied that the well known six- vertex model itself might exhibit such an exotic continuum limit if properly staggered. 
The purpose of this paper is to discuss these results further, and put them on considerably firmer grounds. 

Indeed, the evidence for a continuous spectrum of critical exponents is not so easy to obtain from studies of a finite 
lattice model. What was really established so far in pj [9] was that low energy levels appeared with extremely high 
degeneracies in the limit of long chains, and that naive calculations of finite size corrections indicated truly infinite 
degeneracies. This was — thanks to complementary arguments, such as mappings onto sigma models, or abstract 
construction of WZW theories on supergroups — interpreted as strong indications for a continuous spectrum in the 
scaling limit. Direct evidence was however missing, together with estimates of the measure of integration on the 
continuous spectrum, if any. These issues will be resolved here. 

The paper is organized as follows. In section 2, the critical antiferromagnetic Potts model and the related staggered 
six-vertex model are defined. Symmetries and limiting cases are studied. It is shown in particular that the model 
coincides with the OSP(2\2) model of .7] in one limit, and the 5*0(4) model in another. In section 3, the solution of 
the model using the Bethe ansatz is discussed. In section 4, a detailed analysis of the spectrum of conformal exponents 
from the Bethe equations is carried out. Very accurate evidence for the existence of a continuous spectrum is obtained, 
together with some information on the measure. This information is used in section 5 to relate the results to theoretical 
expectations, in particular those of the supersphere sigma model of [TUj- Elements for a physical interpretation of the 
continuous spectrum are proposed in section 6. Conclusions are gathered in section 7. 

2 The staggered six-vertex model and its integrable structure 

2.1 The general integrable six- vertex model 

On the square lattice C of 2N vertical lines and M horizontal lines, associate the complex number v(J) to the J-th 
vertical line, and h{I) to the /-th horizontal line (see figure[IJ. The parameters v( J) and h(I) are called line rapidities. 
On this lattice, define the general inhomogeneous six-vertex model with local weights given in terms of the difference 
ujj — h(I) — v(J). The weights that satisfy the Yang-Baxter equations are obtained by taking equations (1.5)-(1.6) 
of ref. [llj and performing the substitution : 




(3) 
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Figure 1: The spectral parameter uu denned by the line rapidities. 
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Figure 2: The allowed configurations of the six-vertex model. 



The parameters A/j must satisfy the additional condition : 

A/,j+iA/ + i,j = A^jAz+i^j+i (4) 

The parameters Xjj do not alter the eigenvalues of the transfer matrix, thus they play the role of gauge parameters. 
The parameter A has the value : 

A = — cos 7 (5) 

2.2 The staggered six-vertex model associated to the critical antiferromagnetic Potts 
model 

The anisotropic Potts model on the square lattice C is defined by the partition function : 

^Potts = ^2 J! cx P[Ji^( s i, s j)} )~l exp[J 2 5(si,Sj)] (6) 

{s^ (ij)even (ij)odd 

where each spin Sj lives on a vertex of £ (white circles), and each coupling factor is associated to an edge of C (dotted 
lines). The spin can take Q distinct values. The sum is over all spin configurations, and each product corresponds 
to one type of edge of C. 
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Figure 3: The square lattice £ (clotted lines and white circles) and its medial lattice £ (full lines and black circles). 



When the couplings J\ and J 2 are negative, the model is antifcrromagnetic. In this domain, the critical line 
separating the paramagnetic and antiferromagnetic phases is given by the condition : 



(e Jl + 1) (e Ja + 1) = 4 - Q (7) 

The model defined by eq. ((6|) on the square lattice C can be mapped to a six- vertex model on the square lattice £' . 
This lattice is represented in full lines and black circles, and is called the medial lattice of C (see figure[3]). The weights 
of this six- vertex model depend on the parameters J±, J 2 and Q. We use the notations : 

^/Q = 2cosj (8) 

e Jl - 1 e j2 - 1 



The equivalent six- vertex model has weights a;i, . . . , u)q on the even vertices and wJ,...,Ug on the odd vertices, where : 

ui,...,u e = l,l,ai,ari,e i7 / 2 + a:ie- <T / a 1 e-*T/ 2 +x 1 e i T/ a (10) 
wi,...,o/ 6 = x 2 ,a; 2 ,l,l,e- l ' y / 2 + a ; 2 e^/ 2 ,e^/ 2 +x 2 e-^ 2 (11) 

The parameter A is independent of the parameters x± and a; 2 : 

Wl^ 2 + - W 5 W 6 + ~ WgWg 

A= = ^— ^ — - = — cos 7 (12) 

The criticality condition can be parametrized by : 

s inM cos( 7 - u) 

x\ = r , x 2 = (13) 

sm(7 — u) cosu 

When this condition is satisfied, the weights (p~0|) — (fTTj) of the Potts model correspond to a particular case of the 
integrable six- vertex model ([3]). The rapidity v(J) is equal to (resp. 7r/2) when J is even (resp. odd). The rapidity 
h(I) is equal to u (resp. u + 7r/2) when / is even (resp. odd). This configuration of the line rapidities is shown in 
figure [H The gauge parameter is set to A = e %1 l 2 at every vertex. Thus, the partition function of the Potts model 
at the antiferromagnetic critical point is described by the two-row transfer matrix of the integrable six-vertex model 
with the above choice of the rapidities. We call this matrix T(u) (see figure [5]). 



4 











































/ 


\ / 


\ / 


\ / 


\ 



71/2 7C/2 

Figure 4: The line rapidities corresponding to the antifcrromagnetic critical point of the Potts model 
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Figure 5: The transfer matrix T(u). The value next to each vertex is the spectral parameter uij. 
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Figure 6: The R matrix of the six-vertex model. 
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Figure 7: The Yang-Baxter relation for the six-vertex R-matrix. 



2.3 The 7^-matrix and the thirty-eight-vertex model 
2.3.1 Building block : the i?-matrix of the six-vertex model 

Definition. The i?-matrix of the integrable six-vertex model is defined by its matrix element equal to the 
Boltzmann weight of the configuration shown in fig. [6] Let V be the Hilbert space generated by the vectors | f), | |). 
Then the i?-matrix is a linear operator mapping the space ® V a onto Vp (S> V v . The matrix elements of the R- 
matrix with spectral parameter ujj = u and gauge parameter Ajj = A are the Boltzmann weights ([3]). In the basis 
(|TT),|U),|IT),|U)),thei?-matrixis: 



R(u,X) = 



/ 


sin(7 — u) 


















A e~ lu sin 7 


sin it 












sinu 


A -1 e m sin7 







V 











sin (7 — 


u) j 



(14) 



Symmetries. The i?-matrix satisfies the Yang-Baxter equations shown in fig. [7| and the inversion relation : 

R(u, X)R(~u, A -1 ) = sin( 7 - u) sin(7 + u) 1 (15) 

It also has the symmetry property : 

R(u + 7T, A) = —R(u, A) (16) 



G 



The i?-matrix preserves the total magnetization : 

< + < = °% + vl (17) 



Relation to the Temperley-Lieb algebra. When the gauge is set to A = 1, the i?-matrix can be written in terms 
of the Temperley-Lieb generator E : 

R(u, l) = sin(7-u) i + sinuE (18) 

where : 

/ \ 

e"^ 1 
1 

\ j 

ViNt define the operators 



E 



(19) 



In the Hubert space of row configurations V\ (8 ■ 

E m = tx O • • • <gs l m _i ®E® ll m+2 <8> ■ ■ ■ ® I2JV , m=1...2N (20) 
with a non-trivial action on the space V m ® V m +i- This family of operators satisfy the Temperley-Lieb algebra : 



Q E„ 



Ejjji E nl E nl j r \E m — E m E rn —-iE m 
E m E m > = E m i E m if |m - m'\ > 1 

The operators E m can be expressed in terms of the Pauli matrices : 

1 



E„ 



or, in a more compact form 



S m = (7+ (T m+1 + £7 m cr+ +1 + -(1 - CT^O-^ +1 )e ll < 



(21) 
(22) 
(23) 



(24) 



(25) 



2.3.2 The 7^-matrix. Conservation laws. 



The six-vertex model defined in section [2721 is not homogeneous, and the transfer matrix T(u) is built using i?-matrices 
with different values of uij. One can construct a homogeneous model by considering the 7?.-matrix, acting on double- 
edges. The double-edges live in the Hilbert space : 



V = V ® V 



(26) 



The 1Z- matrix acts on the product space V <£> V. 

As a consequence of the magnetization conservation by R, the 7?.-matrix also preserves the total magnetization 



(27) 



'02 ' ""l 1 ""2 

When the gauge parameter is set to A = 1, another conserved quantity can be constructed. Start from the operator c : 

c = -(cos 7 )- 1 i?(7r/2, 1) = (cos7)- 1 J R(-7r/2, 1) (28) 
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Figure 8: The matrix 1Z(u), denning the antiferromagnetic critical point of the Potts model. 
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Figure 9: A graphical representation of the identity : (c® c) _1 7?.(c<g) c) = 1Z. 



This operator can be expressed in terms of the Temperley-Lieb generator 

c = 11 - (cos7) _1 £; 
According to the inversion relation (|15[) . this operator has the property : 



c z = t 



The 7?.-matrix obeys the conservation rule 



[7£(u),c®c] = 

This is a consequence of the Yang-Baxter equations and the inversion relation, as shown in fig. [5] 



(29) 

(30) 
(31) 



s 



The eigenvectors of the operator c are : 

|0) = (2cos 7 )- 1 /2 ( e <7/2| U ) _ e -n/2| ^ (32) 

|0) =(2cos 7 )- 1 / 2 ( e -^ 2 |U)+e^ 2 ||T)) (33) 

|+) =|TT> (34) 

h) =IU) (35) 

The eigenspace associated to the eigenvalue 1 is {|+), |— ), |0)}, and the eigenspace associated to the eigenvalue —1 is 
{|0». 

2.3.3 Mapping to the 38-vertex model 

The coefficients of the 7?.-matrix in the basis : 

(!+),!-), |0),|0))®(|+) ) |-),|0),|0)) (36) 

define the Boltzmann weights of a 38-vertex model on the square lattice (see figure [TU]). In this vertex model, each 

edge carries an arrow or a thick line. The state |+) {resp. |— )) is represented by an up or right {resp. down or left) 
arrow. The state |0) is represented by an empty edge, and the state |0) by a thick line. Setting : 

7o = 7r - 27 (37) 

u = -2u (38) 

a ,6 ,c = sin(7o - wo),sinu ,sin7o (39) 
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weights of the 38- vertex model read : 
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Figure 10: The 38- vertex model 
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2.3.4 The limit 7 -> tt/2 

If the parameter 7 is set to n/2 and the spectral parameter u is fixed, then the weights of the 38- vertex model are 
equal to zero, except : 

(1) (8) (6) (7) 1 . 2 o 

a\' = a\' = a\' = a\' = — sm 2u 



_ ., _ 1 

t 10 = u lx = u, 12 = u 13 — 



(1,2) _ „d,2) _ „ ( i,2) _ „(u) = : . n22ll 



4 

a 14 = &i5 = ai6 = &17 = — — sin 2 2u 
The 7?.-matrix is proportional to the "graded permutation" P : 

U = - sin 2 2u P (7 -> vr/2, m fixed) (40) 

P|j)|fc) = (-i) Wff *|fc>tf> , i,fee{0,Q,+,-} (41) 
5o,5o,5+,g- = 0,0,1,1 (42) 

This trivial limit can be avoided by scaling the spectral parameter as : 

7T 7T 

7 = - + e , n=- + ew (43) 

where e — > 0~ and u; is fixed. In this rescaled limit, the Boltzmann weights are proportional to e 2 . 

Note that the states |0), |0) become degenerate, but the following combinations remain non-degenerate : 

|0)+t|5) = e ^ 4 [tan(- e /2)]- 1/2 (|U)+ i |4T» (44) 
|0> — i|D> = eT™^ [tan(-e/2)] 1/2 (| \{) - %\ IT)) (45) 

Denote af^ 1 the matrix elements of TZ = 1Z/(—e 2 ) in the basis (|0), |0), |+), |— )). One gets : 

aW =g (8) = ( 1+w _ v f) 

-(2) -(4) 

a\ = a\ = w 

-(3) ~(5) -(1,2) -(1,2) -(1,2) -(1.2) /, x 

a i a i a 2 a 3 — a\ ' = a 5 = (I — w) 

~(6) -(7) -(1,2) -(1,2) -(1,2) -(1,2) v 

a\ = a\ = cl\q = a\i = a-^ 2 = a 13 = w(l — wj 

-(1,2) -(1,2) 

a e — a 8 ' = iw 

-(1,2) -(1,2) 
fl 7 fl 9 

S14 = O15 = (1 - w) 2 
016 = Ol7 = W 2 

5i8 = S19 = (1 - 2io) 

These weights can be related to an integrable loop model with OSP(2\2) symmetry. Indeed, the matrix TZ can be 
expressed as a combination of the identity, the permutation operator P defined in equation (|4ip . and the Temperley- 
Lieb operator E. The latter is defined in tensor notation as a contraction of the spaces V M ,V a and Vp,V v (see 
figure © : 

Kl = -M-^W (46) 
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Figure 11: The three allowed vertices of the dense intersecting loop model. 



where J is the bilinear form in the basis (|0), |0), |+), |— )) 
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tfj = jji = i 



(47) 



By construction, the operator obeys the Temperley-Lieb algebraic relations (j2T|) 

E 2 = Tr(JV) E = 

In the block {|0) ® |0), |0) ® |0), |+) ® |-), |-) ® |+)}, the matrix of the operator E is : 

\ 



H23J), with 



(48) 



/ 1 1 -i i 

1 1 -i i 

-i -i -1 1 

\ i i 1 -1 / 



The matrix 1Z is equal to 



72. = (1 - w)l+ wE + w(l - w)P 



(49) 



(50) 



These weights define an integrable loop model [12], with three allowed vertices, and a loop fugacity q = 0. The 
graphical correspondence is shown in figure [TTJ Note that according to equation (fT5| , the Potts model is isotropic 
when u — -f/2 + tt/4, which is equivalent to w = 1/2. This is also the isotropic point of the loop model (fST)]) . 

We conclude that, in the limit 7 — > 7r/2, the staggered six-vertex model coincides with the OSP(2\2) integrable 
model. In this limit, the arrows represent the fermionic coordinates, and the thin and thick lines the bosonic ones. 
The equivalence in the boundary conditions has to be treated with some care however: the periodic vertex model 
corresponds to the antiperiodic sector of the OSP system, with an effective central charge c c g = l+(— 2) — 24x (— 1/8) = 
2. This means that the OSP symmetry is broken in the periodic vertex model, where the fermions are twisted - i.e. 
they become a complex (Dirac) fermion instead of symplectic fermions. 

2.3.5 The limit 7^0 

The approach for this limit is similar to the case 7 — > n/2. Consider the following limit : 



7 



-we 



(51) 
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where e — ► + and w is fixed. It is convenient to perform the change of basis : 

|0) =t|0) 
|3) =i\o) 



(52) 

(53) 
(54) 
(55) 



Define a transformation of the TvL-matrix that does not affect the partition function : the Boltzmann weights of the 
vertices are multiplied by (—1) for each 7r/2-turn of the thick lines (the weight a\ is not affected). Since the number of 
such turns is even for every thick polygon on the lattice, the partition function is invariant under this transformation. 

Denote a| the matrix elements of the (rescaled) resulting matrix TZ = lZ/(—e 2 ) in the basis (|0), |+), |— )). 
One gets : 



-(1) -(8) 



(1 



s( 2 ) _ s( 4 ) 



.(3) -(5) -(1,2) -(1,2) -(1,2) -(1,2) 
a l a l a 2 a 3 a 4 a 5 

= (6) _ =( 7 ) - sf 1 ' 2 ) - ^(L 2 ) _ = (1,2) _ -(1,2) _ 
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The corresponding loop model is constructed with two generators P and E. The permutation operator P is defined 
as : 

Pfivtia = 5p a 5 vtl (56) 
The Temperley-Lieb operator E with parameter i/Q = 4 is built using the simple contraction : 



E. 



J i>(3 J a 



where 





( 1 





°\ 


J = 





1 











1 
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o / 


K = 


(14 


-to)l- 


wE 



J =1 



(57) 



(58) 



The matrix 1Z is equal to 

+ w(l + w)P (59) 

These are the integrable weights of the loop model defined in [T^] with a loop fugacity q = 4. Note that, in this 
regime, the Potts model is anisotropic for any value of to. The loop model itself is isotropic when w = —1/2 : the 
relative weight for loop crossings is then equal to —1/2. Following [T2], this model is a graphical version of the SO (A) 
integrable model based on the vector representation. Using that so(4) = si (2) + sl(2) one can expect it to decouple 
into two copies of the isotropic six-vertex model or XXX spin chain, a feature we will confirm when discussing the 
associated hamiltonian or Bethe equations. 
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2.4 Transfer matrices 

In this section, we set the gauge to A = 1, so that the transfer matrices can be related to the generators of the 
Temperley-Lieb algebra. 

One-row transfer matrix. The two-row transfer matrix T(u) is the product of two one-row transfer matrices : 
T{u) = T(u + n/2) T{u). (see fig. [5J . As a consequence of the Yang-Baxter equations and the inversion relation, when 
periodic boundar conditions are imposed in the horizontal direction, the one-row transfer matrices commute : 

[T(«),T(«')]=0 (60) 

The one-row transfer matrices have the property : 

T(u + tt/2) = (-l) N e ip T(u)e- iP (61) 

Define the operator e~ %p as the shift of all vertical edges by one site to the right. The operator c defined by equation (f28|) 
is used to build a conserved quantity of the transfer matrix. Denoting c m the charge operator acting on the space 
V m ® V m+ i, the global charge is : 

C = ci x c 3 x • ■ ■ x c 2N -i (62) 
C 2 = 1 (63) 
If u = 0, the transfer matrices T(u), T(u + tt/2) become : 

n \ N 

T(0) = f -sin2 7 J e~ ip C (64) 

/ 1 \ N 

T(tt/2) = I --sm2 7 J C e~ lP (65) 
The matrices (|62p . (|55|) are represented in figure [T2l 

Conservation laws for T. As a consequence of the conservation of the total magnetization and the "charge" c by 
the 7?.-matrix, the transfer matrix T preserves the total magnetization : 

S=\k + --- + ^n) (66) 

and the "charge" C defined above. Since 7" (it) also commutes with the transfer matrix T(0), T(u) is invariant under 
the action of the one-site shift operator e~ lP . 

2.5 Hamiltonian limit of the transfer matrix 

To compute the derivative of the matrix T(u) at the points u = 0,u = n/2, it is convenient to write this matrix in 
terms of the i?-matrix : 

T 0i...faN( u \ = y R^( U )R^( U - W /2)...R^' 3 ^-] (u)R^™ (u-n/2) 

Denote by 5 the derivative with respect to u. According to eq. Q18[l : 

5R(u) = — cos(7 — u)l+ cos u E (67) 
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(c) 

Figure 12: Diagrams defining the transfer matrices C (a), e~ tP C (b) and Ce~ lP (c), for N = 4. The white circles 
represent c operators. 



Differentiating the previous expression yields : 

AT-l AT 

sin 5t 

2 



/l \JV-1 JV 

«X0) = r sin 2 7 cos 2 7 i 2j _i + sin 2 7 %) (68) 



where the matrices t 2 j-i and <2j ar e defined by the diagrams on fig. 1131 Using the graphical representation of the 
matrices T(0), t 2 j-\, faj in figures [T2l TT3l and the property ([30)) . one gets the intermediate results : 

[T(0)]~ l t 2 j-i = (isin27) %_ 3 «%-! ^j^TOj-! %- 3 ®^-2 ( 69 ) 



[^(0)] % - Q sin2 ^) cv^-i®^ (70) 
The Hamiltonians i?2 are defined as the logarithmic derivatives of the transfer matrix T{u) at the points u = 0, 7r/2. 

Hx = i sin 27 [r(o)]-Mr(o) 

TV 

= cos 2 7 11 + cos 7 + E 2j ) - {E 2j E 2:j -i + E 2 j-iE 2j )] 

3 = 1 

(71) 
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Figure 13: Diagrams defining t 2 j-i and t 2 j of eq. (|68p. for TV = 4 and j = 3. The white circles represent c operators. 
The arrow points to the vertex that differs from T(Q). 



Using eq. ([61 



H 2 = -sin2 7 [T(Tr/2)]- 1 ST(ir/2) 

iP tt -iP 

= e Hie 

N 

= ^[-cos2 7 ]l + cos 7 {E 2j + E 2j+1 ) - {E 2j+1 E 2j + E 2j E 2j+1 )] 

3 = 1 



The Hamiltonian associated with the two-row transfer matrix T{u) is defined as : 

W=ism2 7 [T(0)]-MT(0) 

According to the commutation relations (|60p . 

U = H x + H 2 

2N 

= ^ [- cos2 7 1 + 2 cos 7 E m - {E m+1 E m + E m E m+1 )\ 

m— 1 

Note that the Hamiltonian Ti is the sum of two commuting parts : 

H = H 1 +H 2 , [H 1 ,H 2 }=0 
A consistent decomposition of the momentum operator is : 

e- 2lP = (e- lP C) x {Ce- lP ) = (C e - lP ) x ( e -' iP C) 
Indeed, the operators Hi, H 2 , (e~ lP C), (Ce~ lP ) all commute with one another. 



(72) 



(73) 



(74) 



(75) 



(76) 
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Using the expression (|25[) of the Temperley-Lieb generators, the Hamiltonian ([74"]) is expressed in terms of Pauli 
matrices : 



n 



2N ^ 
^2 [ ~~ ( a m a m+2 + a m, (T m+2) + sin2 7 ^m^m+l ~ ^^m^ 



m+2 



m u m+1 / 



• cos 2 7 1] 



(77) 



In the limit 7 — > 0, the Hamiltonian (|77|) describes, as expected from the 50(4) identification, two decoupled 
ferromagnetic XXX spin-chains : 



1 N 1 N 

H = --Y^2j-l ■ ^2j+l + 1) - - • ^2 i+ 2 + 11) (7 = 0) 



(78) 



3 Bethe equations, Bethe states and eigenvalues 

When periodic conditions are imposed in the horizontal direction, the staggered six-vertex model is solvable by Bethe 
ansatz [11] . We call r the total number of particles. 

3.1 Bethe equations for two types of particles 

The Bethe ansatz equations are : 



Mj S {1, . . . , r} exp [2iNk((Xj)] = - exp [-i<p{a } , a{)] 



The one-particle momentum and the scattering amplitude are given by : 

sinh(a + 17) 



exp [2i/c(a)] 
exp [i<t>(a, a')] = 



sinh(a — ij) 

sinh h[a — a' — 2i 7 ) 
sinh h(a — a' + 2ij) 



(79) 



(80) 



(81) 



The roots ctj describing the ground state and the physical excitations are expected to sit on the two lines Im(a) = ±7r/2. 
Define two types of particles : 



a 



Aj+wr/2, j = l, ...r+ 
Hj-iw/2, j = l,...r- 



The one-particle momenta are 



exp 



2ik(X) 



exp 



2ik( \±i- 



The scattering amplitude between two particles of the same type is : 



r / 7T 7T 

exp [i<f>i(\, A')] = exp icj> \\ ± i — , A' ± i— 



cosh(A + 17) 
cosh(A — iy) 



sinh i (A - A' - 2i-y) 
sinh i (A - A' + 2i 7 ) 



(82) 
(83) 

(84) 
(85) 
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The scattering amplitude between two particles of different types is 



r / 7T 7T 

exp [itp^i(X, A')] = exp i<p y\ ± i— , A' =p i — 



cosh i (A - A' - Hi) 
cosh | (A - A' + 217) 



Define the shifted scattering amplitudes as odd functions of (A — A') : 

exp [i6i(A, A')] = - exp [%(A, A')] 

exp [j©_i(A, A')] = exp [i0_x(A, A')] 
The Bethe equations (|79|) split into two sets : 



(86) 



(87) 



-i©_i(Aj 



1=1 



;=i 



Vj S {!,..., r_} e 2iNk ^ = {-l) r — 1 JJ e -^-ife,A ( ) 



iOj A) tt -»ei(/i3,Mi) 



Take the logarithm of equations (|59")) ~ (f9"0"f : 

Vj e {1, . . . ,r+} 2M(A,) = 27T/+ - ^e^A,, A ; ) - E Q^{\ 3 , W ) 



i=i 



z=i 



Vj € {1, . . . , r_} 2Nk( f x j ) = 2ttIt - £ Aj) - E Gj^, W ) 



(89) 
(90) 

(91) 
(92) 



i=i 



i=i 



The "Bethe integers" 7? follow the following rules : if r + (resp. r_) is even, then all the Tt" (resp. Ij ) are half-odd 
integers; if r + {resp. r_) is odd, then all the 1^ (resp. Ij) arc integers. Summing equations (|9ip ~ (|92[) over j and 
recalling that the functions 0±i are odd, one relates the total momentum to the Bethe integers : 



fctot = £ fc(A,) + E fc(Mi) = ^ (E^+E 7 : 



i=i 



The one-particle momenta and the scattering amplitudes can be written : 



tan 



fc(A) =tanhA tan 7 



tan 



tan 



9i(A,A') 



e_i(A,A') 



tanh - (A — A') cotan7 



tanh-(A — A') tan 7 



(93) 

(94) 
(95) 
(96) 
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3.2 Bethe states 

Form of the Bethe states. The one-particle states ip a are "inhomogeneous plane- waves" , represented in real space as : 

(p a (2m-l) = coshi(a-i 7 ) e 2lfc ( Q ) m (97) 
ip a (2m) = - sinh 1 (a + i 7 ) e 2lk ^ m (98) 
with m = 1, . . . , N defines the position of the particle. The general Bethe state is given by the linear combination : 

<p ai ,...,a r (Zl, ■ ■ ■ , X r ) = ^2 A Hu...,Prfa P1 (xi) ■ ■ ■ fa Pr (x r ) (99) 

p 

where the sum is over all permutations P — {pi,P2> ■ ■ ■ ,Pr} of the integers 1, . . . ,r, and x\, . . . , x r are the positions of 
the particles. The coefficients A plt ___ tPr are given in terms of the scattering amplitudes Sjj : 

Api,...,p r = e pY\_ s Pj,%>i ( 10 °) 

i<j 

Sjj = sinh \{cxj — a; — 2ij) (101) 

where ep is the signature of the permutation {p\, . . . ,p r }- In particular, when the roots ctj sit on the two lines 
Im(a) = ±7r/2, we introduce another notation for the Bethe state : 

y(Ai,...,A r+ |jti,...,/*r_) = f\ 1 +m/2,...,X r+ +m/2,n 1 -m/2,...,iJ, r _-i-ir/2 (102) 



Symmetries. By construction, the states (p ai ,...,a r are antisymmetric under the exchange of the otj. A shift ctj — > 
Q.j + 2i7r results in a phase factor : 

(103) 



Action of C . The action of the operator C on the Bethe states is a shift ctj — > ay + i7r : 

C(p ai ,..., ar = {-i) T <Pai+i-K,...,a r +m (104) 

The operator C, acting on the states (|102j) . exchanges the two lines : 

C<P(x u ...,x r+ \» u ...,» r _) = e^ {r ^ r - ] 5P0' 1 ,..., /lr _|Ai,...,A r+ ) (105) 
A special class of states (to be discussed in section 13. 5|) is defined by the following constraint on the Bethe integers : 

r + = r_ = r (106) 

Vje{l,...,r} , // /, /, (107) 
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These states are eigenvectors of the operator C, with eigenvalue one : 

CV(Ai,...,A r |Ai,...,A r ) = ¥>(Ai,...,A r |Ai,...,A I .) 

The other eigenvectors of C are given by the linear combinations : 



l,...,A r+ |,il,-,A»r_) ~~ ^ ( ¥ '( Al >-' A '-+l' tl '->' i '--) ± 6 : 



■(r+-r_) 



i,...,/V_ |Ai,...,A r+ )^ 



^(Ai,...,A r+ |jii,...,j*r_ ) ~' = ^'(Ai,...,A r+ 



Action of the shift operator. The shift operator e iP has a similar action on the Bethe states : 



— 7 P 

e </>ai, ...,o 



J^J ^ e -*pi(«i+* 71 ') 



where : 



exp [ipi(a) 



sinh 2 (a + 17) 
sinh 5(0! — 17) 



3.3 Eigenvalues 

The eigenvalue of the transfer matrix T(u) associated to the Bethe state (p ai ,...,a r isi '■ 

A(ai, . . . , a r \u) = fi(ai, . . . , a r \u)ii(ai + in, . . . ,a r + m\u) 

where 

/i(ai, . . . ,a r |u) = (i/ 2 )^ x 



|sin2(u - -f)] N f] sinh ^( Q j + (- sm 2 ?i ) w TT sinh 5( a J _ 2 * u + 3i ^) 



l sinh 2 (ctj — 2iu + 27) 



^ sinh 5 ("j ^ 2iu + i'y) 



In the limit u — > 0, the energy of the state tp ai ,...,a r is defined as : 

cv \ 1 • O ^logA i n x 

t(ai,...,a r ) = -sin 27 — (ai, . . . ,a r |0) 



= -2ATcos27 + V' 



2 sin 2 27 



cosh 2a, — cos 27 

By construction, the state ip ai ,...,a r is an eigenvector of the Hamiltonian H, with eigenvalue -E(ai, . 
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3.4 Reminder : the homogeneous six-vertex model 

Consider the homogeneous six- vertex model defined by the Rq(uq) matrix given in equation (|14p , with parameter 70, 
on a lattice of width N. 

In the anisotropic limit Uq — > 0, the logarithmic derivative of the transfer matrix T (u ) is equal to the XXZ 
Hamiltonian : 

H = -sin 7o [T (0)]- 1 ^(0) (117) 

du 

N 

= ^ (cos 7o -O (118) 

771 = 1 

where the operator is a generator of the Temperley-Licb algebra with parameter 70 , like in equation (|24p : 

Ki = \ K °m+i + < vfn+i - cos 7o«77 <+i - I) - * sin 7o « - a z m+1 )] (119) 

On a lattice with periodic boundary conditions, the system is solvable by Bethe ansatz. The Bethe equations for 
r particles are : 

r 

Vje{l,...r} exp[iiVfc (A,)] = -[]exp[-i0 o (A,,A i )] (120) 

i=i 

where the one-particle momentum and the scattering amplitude are : 

e X p[i fc0 (A)] = Si ^! 70 :^ (121) 
sinh(|7o + A) 

r , . /m sinh(A — A' + i7n) 
exp [i<po(\, A' = . J, . 122 

smh(A — A' — i-fo) 

The associated eigenvalue of the transfer matrix T (uo) is : 



., , r • 1 r , smh(A 3 + iuq + ^70) 

A (Ai,...,A r \uo) = sm(7o-M ) I I — r-TTT — ~ 1 — T 

XJL » smh(A j + iu - 570) 



( ■ \ N TT ( smh ( A i + m o - TTo) 
(srnuo) ttt; : 1 — T 

y smh(Aj + iu - 270) 



(123) 



3.5 Common eigenvalues between the staggered and homogeneous six-vertex models 

If the parameters of the staggered and the homogeneous six- vertex models are related by : 

70 = tt- 2 7 (124) 
m = -2m (125) 



then one has the relations 



exp 



2ik{\) = exp[ik (\)) (126) 
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fc(A,A')+^-i(A,A') = o (A,A') (127) 

These relations suggest that the states : 

V(A ll ..,A,|Ai,...,A r ) (128) 

have properties described by the homogeneous six- vertex model. The states (|128|) are obtained as the solution of the 
Bethe equations (f9Tj) -(f92 |) when the Bethe integers on the two lines are the same : 



r+ =r- =r (129) 

VjG{l,...,r} , /+ = /-=/, (130) 

As a consequence of relations (|126l) - (|127[) . the Bethe equations (|9T|) - (|9^1) are equivalent to the Bethe equation (|120p 
for the homogeneous six-vertex model. Moreover, the eigenvalue of the transfer matrix T(u) can be written in terms 
of the eigenvalue (|123| : 

A (Ai + i^, Ai - . . . , A r + Ar - i||u) = (-1/4)* [A (Xi, K\u )} 2 (131) 
As a consequence, in the anisotropic limit, the energies £, £ °f the Hamiltonians 7i, Hq are related by : 

/ 7T 7T 7T 7T\ 

£ V 1 + * 2 ' Al ~~ * 2 ' ' ' ' ' Ar + i 2 ' Ar ~~ * 2 J = 2f °( Al ' ' ■ ' ' Xr > ( 132 ) 
The total momenta are equal : 

hot = fco.tot (133) 

4 Finite-size study of the Bethe equations 

4.1 Reminder : the XXZ spin-chain Hamiltonian. 

4.1.1 Definition 

The XXZ Hamiltonian (|118|) on a lattice of width N can be written : 

N 

H = -1/2 + «+i ~ cos 7o + 1)] (134) 

m— 1 

< 70 < vr (135) 

with periodic boundary conditions : 

<+i = < (136) 

4.1.2 Ground state 

The ground state of the Hamiltonian (|134[) is given by the symmetric "half-filled Fermi sea" in the sector with r = N/2 
particles. Consider one of the particle-hole excitations around the Fermi momenta ±kp ■ 

k F - feo/2 -> k F + k /2 (137) 

— kp + ko/2 — > —kp — fco/2 (138) 
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These excited states have total momentum ±ko. In the thermodynamic limit N — ► oo with the ratio r/N fixed, the 
Bethe equations (|120p give a linear integral equation for the density of particles. Solving this equation gives access to 
the dispersion relation of the excitations (|137p — (|138[) : 

£ (k ) ~ v \k \, v = n Sin 70 (139) 

7o 

This shows that, in the thermodynamic limit, the spectrum has no finite gap above the ground state. This is an 
indication that the equivalent (l+l)-dimensional quantum field theory is conformally invariant. The velocity vq 
appears in the finite-size determination of the central charge : 

£ 0tSr (N) = N eo.oo - 9^ + < N '^ ( 14 °) 

The central charge of the Hamiltonian Hq is c = 1. The energy density per site is : 

2 f°° dx 

e ,oo =cos7o-2sm 70 / — — — r- — -— r (141) 

J Q cosh(7rx)[cosh(27ox) — cos 70 J 

4.1.3 Low-lying excitations 

A class of Bethe states are low-lying excitations, with energies : 

£o(N) = W^V) + + oiN- 1 ) (142) 

where x is called the physical exponent. We describe the primary states, and then their descendants. In the sector with 
r = N/2 — n particles, the lowest-energy state is given by the distribution of the Bethe integers which is symmetric 
around zero. In the same sector, denote by {n, m} the state obtained after m backscatterings from the left Fermi level 
to the right Fermi level : 

h,...I r = - r -^-+m 1 .. r -^-+m (143) 
The physical exponents of this class of states are : 

4°L=« 2 f+- 2 f, 9= n —^ (144) 
2 2g 7r 

Descendant states are obtained by performing particle-hole excitations on the above states. For example, starting 
from the ground state distribution and setting I r to (r + l)/2, one obtains a state with physical exponent x = 1. 

4.2 Ground state of the Hamiltonian H 

In this section, the ground state of the Hamiltonian (|74|) . corresponding to the anisotropic limit of the staggered 
six- vertex model, is discussed. The system width N is assumed to be even. The case of an odd system width N will be 
discussed in section 14.41 The ground state of the Hamiltonian H is given by the symmetric distribution of the Bethe 
integers in the sector r + = r_ = N/2 : 

I+ I+ _N/2-l_N/2-l N/2-l 

N/2 - 1 N/2 - 1 N/2 - 1 . . 
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Note that this state has twice the energy of the ground state of Hxxz ■ 

Consider the double particle-hole excitation obtained by performing the transformation (|137[) on both Fermi seas. 
This excitation is of the type (|129p - (1130[) . and therefore, using equations (|139p . (|133|) and (|132p its momentum and 
energy are : 

k = k (147) 
£{k) ~ 2 vo \ko\ (148) 
Thus, the rapidity of the particle-hole excitations around the ground state of the Hamiltonian Ti. is : 

!j = 2u (149) 

Compare the asymptotic behaviors of the finite-size ground state energies for Hamiltonians H , Ti : 

£o, gr (N) = Ne 0<oo - |^ + o{N- 2 ) (150) 

S gI (2N) = 2Ne 00 --^^ + o(N- 2 ) (151) 

As a consequence of the identities (|132|) and (I148p . the central charge of Hamiltonian Ti is c = 2. The energy density 
per site of Hamiltonian Ti is : 

eoo = e ,oo (152) 

4.3 Low-energy spectrum of the Hamiltonian TC 
4.3.1 Bethe excitations 

The excitations considered are combinations of particle and backscattering excitations on the Bethe integers Ij ,IJ . 
Denote tp( n+ ,n-),(m + ,m-) the Bethe state defined by the numbers of particles : 

r+ = N/2-n+ (153) 
r_ = iV/2-ri_ (154) 



and the Bethe integers 



r + — 1 r_i_ — 1 

IJ + = h m+, . . . , — - — -I- m+ (155) 



_ _ r_ — 1 r_ — 1 
I l ,...I r _ = hm_,..., — hm_ (156) 



The states described in section 13751 are of the type ^(n.nJXm,™)- Their physical exponents are : 

X(n,n),{m,rn) = 2 X^ m = n? Q + m 2 - (157) 

9 

where 

9 = 27 A (158) 
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The analytical study [9] of the thermodynamic limit of the Bethe equations suggests the following general form for 
the physical exponents : 



+ ix( 7 ,iV)(n + -n_) 2 + i] ^ : ^(m + -m_) 2 



(159) 



where the coupling constant if (7, N) tends to zero as N goes to infinity. One of the key questions is to determine the 
way this constant actually vanishes. 

4.3.2 General structure of the low-energy spectrum 

Note that the states with m + =/= m_ are not part of the low-energy spectrum. Thus, in the following, the discussion 
will concern only the states with m + = m_ = m. For these states, the total magnetization and the total momentum 
read : 

S = (160) 

7T 

fctot = —m{N-n + -n-) (161) 

Equation (| 1 59[) determines the structure of the low-energy spectrum. The quantities {n + + n_ ) and m define a 
sector of given total spin S and total momentum fc tot ■ We call "floor states" the lowest-energy states of these sectors. 
These are the states ^(„. n ),(m,m) an d ^ +1 n j < mm y wnere n,m are integers. The physical exponents of the floor 
states are finite in the limit N — > 00 : 

2 m 2 

X(n,n),(m,m) = 9 n H ( 162 

g 

X(n+l,nUm,m) = t( 2 «+ 1 ) H 1 : ( 16 3) 

4 g 4 

Note that the states <P( n ,n),(m,m) arc those described in section [3~51 Within each sector, higher energy states are 
obtained, starting from the floor state, by "moving" nf particles from the line Im(a) = ir/2 to the line Im(a) = — tt/2 
(or the reverse) . The physical exponent of the resulting state differs from the floor exponent by a quantity proportional 
to the coupling constant if (7, N) : 

X[n+n',n-n'),(m,m) = ^(n,n),(m,m) + ("') 2 -^(7> N ) ( lg 4) 
V(n+n'+l,n-n'),(m,m) = ^(n+l,n),(m,m) + n> {n' + l)K(j, N) (165) 

Thus, if the constant K indeed vanishes in the limit N — > 00, the gaps between the floor state and the higher states 
in the sector should also vanish. The structure of the spectrum is illustrated in figure 1141 

Numerical calculations are used to confirm the form (|159p of the physical exponents, and to determine the scaling 
law for K(j,N). 

4.3.3 Numerical calculation of the physical exponents 

Numerical procedure. When the Bethe integers 1^ are fixed according to equations (|155[) - (|156|) . the Bethe equa- 
tions (|91D (192p are a set of non-linear equations for the variables Xj , [ij . These equations can be solved numerically, 
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'(2,-1), (0,0) 
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Figure 14: The first levels of the fundamental sector (n + + n_ = 0, m = 0) and the sector defined by (n + + n_ = 
l,m = 0). The floor states (including the ground state) are represented by bold lines. The energies are rescaled as 
physical exponents. 
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using the multidimensional Newton- Raphson method [13]. The following starting point for the algorithm was found 
empirically to lead to a good convergence : 



A; nlt = Arctanh 



tan(7r/+/A0 
U 



(166) 



where the additional parameter is set to U = 10. The exponents are then estimated, relatively to the ground state, or 
to the floor state, using : 

£ a -£b ~ ^^-{Xa-Xb) (168) 

where a, b denote any states of the spectrum. 

Results for the floor exponents. These are shown in figures I15H17I These results agree with the form (I159j) . 

Results for the gaps inside a given sector. These gaps are given in equations (|164p ~ (|165p . as a function of the in- 
teger n' and the coupling constant if (7, A). The first step is to check the dependence on n'. See figures [T8H2T1 

Results for the coupling constant. Equations ()164j) - (|165|) relate the gaps within a given sector to the coupling constant 
A (7, A). These are used to determine the scaling law of this quantity. Theoretical arguments (see section [SJ predict 
the following form : 

K ~ -, ^ (169) 

[B + logAf V ' 

where the exponent p may take the values p = 1,2. The exponent p is estimated numerically, assuming that K behaves 
as equation (|169p (see figure [22]) . A crossover is observed between the behaviours : 

K ^> N) * m< VZV (7<?r/2) (170) 

[Byy) + log A] 

A' 

K(tt/2,N) ~ (171) 

v 1 ' ' B' + logA v ; 

The factor ^(7) in the scaling law (|170|) is estimated numerically by eliminating the term .8(7) between two system 
widths. The finite-size estimators are defined as : 

A(7,iV 1 ,iV 2 )=(^ (7)iNri)] _ 1/2 _ [i , (7)iV2)] _ 1/2 ) (172) 
The numerical results (see figure |2"5|) allow us to conjecture the following form for the factor ^(7) : 

A{l) - (173) 

7T — Z7 

The correction term .8(7) depends on the sector considered. 
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Each sector determines a specific scaling law for K{^), N), and these determinations can be compared. In addition, 
the coupling constant if (7, AT) appears in the expression of "infinite physical exponents" : 



"(n,n),(l,-l) — x (n,n),(0,0) 



(174) 



The results are shown in figure [2"4l 
4.4 The case N odd 

When the system width N is odd, the structure of the spectrum exhibits some differences from the case N even. 
The general formula for the finite-size corrections is identical to the formula for N even : 

TTV 2tTV 

£(. n+ ,n-),(m+,m-)(2N) = 2Ne °o - q x 2N X 2 + lyy X (n+,n-),(m+, m -) + o(-/V _1 ) (175) 

where the exponent xr n n _\/ m m \ is given by equation (|159[) . The subtlety is that equations (|153[) — (|154[1 imply 
that n + ,rt_ are half-odd integers. 

The lowest-energy states of the spectrum (|175[) are <p^y 2 _i/ 2 ) (o ov Note that these states belong to the sector 
5 = 0, and are the most closely packed to the imaginary a axis in this sector. The energy of the ground state is : 

£ gI (2N) = 2N eoc ~ x c + oiN- 1 ) (176) 
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Figure 16: The floor exponent znn n _i) as a function of 7, for a system width N = 512. The expected value is 
=9 + 9~ X (equation (152])). 
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Figure 17: The floor exponent a/(„ +1 n N , \ = :E(n+i,n),(o,o) — £(1,0), (0,0) as a function of 7, for a system width N = 1024. 
The expected value is x', n+1 n x , n = gn(n + 1) (equation (|163|) ). 
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Figure 18: The exponent £( n ,-n),(o,o) as a function of n, for parameter 7 = 0.37T. The expected values are xt, 
K(j,N) n 2 (equation (TMjl ). 
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Figure 19: The exponent ^(i+ n ,i- n ),(o,o) = x (i+n,i-n),(o,o) ~ ^(1,1), (0,0) as a function of n, for parameter 7 = 0.37T. 
The expected values are n),(o,o) = 

if(7,iV) n 2 (equation (fl64"ll ). 
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Figure 20: The exponent £(i+„.i-n).(i,i) = £(i+n,i— n),(i,i) ~ ^(ljlJJi,!) as a function of n, for parameter 7 = 0.4-7T. 
The expected values are £(i+ n ,i-ra),(i,i) = ^(j,N) n 2 (equation (|164[) V 




Figure 21: The exponent S(i+n ! -„).(o,o) = x (i+n,-n).(o,o) ~~ a; (i.o),(o.o) as a function of n, for parameter 7 = 0.37T. The 
expected values are a^i+^-n^o.o) = K(j, N) n(n + 1) (equation (|165p V 



32 



1 - iV = 256 - 1024 — I — 
N= 512 - 2048 — X— 

0.8 I 1 

0.1 0.2 0.3 0.4 0.5 

7/tt 

Figure 22: Estimation of the exponent p in the scaling law (|169[) . The estimator for p is obtained by eliminating the 
unknowns A, B from the equations log if (7, N) — log A — p\og(B + log AT), for three system widths. The estimators 
converge slowly to p = 2 for 7 < tt/2, and to p = 1 for 7 = tt/2. 
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Figure 23: Determination of the factor ^(7) in the scaling law (|170p . using a numerical estimation of exponent 
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Figure 24: Various determinations of the factor ^(7), for system widths N = 256 — 512. 



where the effective central charge is : 

2 = 2- 12x (1 /2 : _i/2) : (o,o) 



£(1/2,-1/2), (0,0) 



(177) 



The effective central charge is estimated numerically using equation (|176|) and the exact expression (|141|) of the energy 
density eoo (see figure l25j) . The scaling law for the coupling constant K(^,N) is consistent with equation (|170|) . It is 
convenient to define the physical exponents of the excited states with respect to the ground state : 



x (n + ,n-),(m + ,m-) = x (n + ,n_ ) , (m+ ,m_ ) — a; (l/2,-l/2),(0,0) 

The physical exponents of the floor states are : 



x (n+l/2,n- 1/2), (m,m) 



x (n-l/2,n-l/2) : {m,m) 



"111 



m 
9 



(2n- l) 2 + 



m 
9 



K 

T 



(178) 

(179) 
(180) 



In particular, the states ¥>(-i/2,-i/2),(o,o)> ^(1/2, 1/2), (0,0) have an energy which is twice that of the degenerate ground 
state of the XXZ spin-chain (|134[) with an odd lattice width [14]. These two states appear as excited states of the 
Hamiltonian 7i. 

The gaps inside a given sector are : 



x (n+n' + 1/2, n-n'- 1/2), (m,m) 
*£(n+n' — 1/2, n— n' — 1/2), (m,m) 



^(n+l/2,n-l/2), (m,m) + Kn' (ri + 1) 
Z(n-l/2,n-l/2), (m,m) + K{n') 2 



(181) 
(182) 
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Figure 25: Determination of the factor A(j) in the scaling law (|170|) by the effective central charge c for N odd. 
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Figure 26: The floor exponents 5(„ +1 / 2 ,n-i/2),(o,o) as functions of 7, for a lattice width iV = 1023. The expected 
values are S(„+i/2, n -i/2),(o,o) = .9™ 2 (equation Q179P ). 
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Figure 28: The exponent Z(i/ 2 +„,-i/2-n),(o,o) = *(i/2+ n ,-i/2-n),(o,o) -£(1/2-1/2), (0,0) a s a fu nction of n, for parameter 
7 = 0.3-7T. The expected values are X(i/2+n.-i/2-n),(o.o) = K(j,N) n(n + 1) (equation (|181I0 . 
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5 Interpretation and relation to non-linear sigma models 

As discussed in section 2, the limit 7 — ► tt/2 of the staggered six- vertex model coincides with a particular point of the 
OSP(2\2) loop model of [7]. This is a good starting point to understand the emergence of a continuous spectrum of 
critical exponents. 

5.1 A reminder on intersecting loop models and Goldstone phases 

It turns out that the Mermin Wagner theorem forbidding the spontaneous breaking of a continuous symmetry in two 
dimensions does not hold for supergroups (because of the lack of unitarity), and that models with orthosymplectic 
OSP(m\2n) symmetry do exhibit a low temperature phase with spontaneous broken symmetry provided m — 2n < 2. 
More precisely, consider the non linear sigma model with target space the supersphere S m ' 2n — OSP(m\2n)/OSP{m — 
l|2n), a "supersymmetric" extension of the usual O(N) sigma model. Use as coordinates a real scalar field : 



and the invariant bilinear form 



4> = (01, ■ ■ ■ ,4>m,i>l, ■ ■ -Aln) 
4> ■ <t>' = ^ §a<\> a + ^2 J aP^a^0 



(183) 
(184) 

^ . The unit supersphere is 
(185) 



where J a p is the symplectic form which we take consisting of diagonal blocks 
defined by the constraint : 

6-6=1 



1 

-1 
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Figure 30: Interaction vertices in the OSP(m\2n) symmetric sigma model on the square lattice. 



The action of the sigma model (conventions are that the Boltzmann weight is e 



s ) reads : 



S = ^- [ d 2 xd^ ■ 



(186) 



'a J 




(187) 



The model for g a positive thus flows to strong coupling for m — 2n > 2. Like in the ordinary sigma models case, 
the symmetry is restored at large length scales, and the field theory is massive. For m — 2n< 2 meanwhile, the model 
flows to weak coupling, and the symmetry is spontaneously broken. One expects this scenario to work for g a small 
enough, and the Goldstonc phase to be separated from a non perturbative strong coupling phase by a critical point. 

It is easy to suggest lattice models whose long distance physics is described by the supersphere sigma models. It 
was shown in [Jj that simply taking a square lattice with the fundamental representation of the OSP(m\2n) on each 
link and Heisenberg coupling at every vertex would suffice. More generality, it is convenient to represent the link states 
by lines carrying a label a = 1, . . . , m or a = l,...,2n and express the interactions in terms of the three invariant 
tensors of the algebra, corresponding to the three diagrams in figure 1301 

A trivial rescaling of the partition function and isotropy of the model allows one to set the first two weights equal 
to 1. The remaining weight is a free bare parameter. It was found in [7] that for any w > and to — 2n < 2 the lattice 
model is indeed described at large distance by the UV limit of the sigma model, ie a set of to — 1 free uncompactified 
bosons and n pairs of symplectic fermions, with a central charge c = to — 2n — 1. The value w = is critical, and was 
argued in [JJ [TS] to correspond to the critical point of the sigma model. Note that for m — 2n fixed, there exists a 
particular value of the crossing weight w where the model is integrable, and coincides with the one in |12j . 

Of course, the geometrical representation of the invariant tensors allows for a full geometrical formulation of the 
model, that can then be extended to values of to, n not integer. The model so obtained is made of loops covering 
every link of the lattice, and possibly self-intersecting once on the vertices, with a fugacity to — 2n, and a weight w per 
crossing. The continuum limit was found to be a naive extrapolation from the results at to, n integer. Note that this 
model differs from the usual formulation of the Q-state critical Potts model with m — 2n= ^/Q only in that crossings 
are allowed. This however changes the universality class deeply. For instance, all the i-leg operators which have non 
trivial, Q-dependent scaling dimensions in the Potts model for —2 < \[Q < 2, now have logarithmic correlators with 
vanishing effective dimension. This is because the symmetry being spontaneously broken the fundamental field has 
non vanishing expectation value, and thus the fields <p a do exist in the conformal field theory, unlike in the case of a 
compact boson. 

[Paths that can self intersect at a vertex but not at a bond are called trails in the literature, and "bond self 
avoiding models" by contrast with "site self avoiding" . The case m — 2n = we consider below would correspond 
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to "fully packed trails". People have considered ordinary (dilute) trails which are in the same universality class as 
dilute SAW and trails at the theta point. Those have been shown to be in the same universality class as "growing 
self avoiding trails" , themselves a particular case of the trajectory of a particle moving on a lattice with random 
distribution of scattering rotators. See [THE]. In the context of such trajectories, it has been remarked already that 
when one dilutes the set of random rotators from maximum concentration Cl — Cr — 1/2 (i.e., add intersections), 
the exponents change.] 

5.2 Coupling constant and physical exponents of the OSP(2\2) model 

We now specialize to the case of interest here, the OSP(2\2) model. The UV limit is easy to understand. We can 
parametrize the supersphere 

(<PlJ \<>2 

by setting 



2 ' ^ ^ 2 ' 2V>iV>2 = 1 (188) 



n — 

h = 



COS0 (1 - ipifo) 

sin0 (1 - Vi^a), (f> = (f> + 2Tr (189) 



The action then reads 



S=^- J d 2 x [{d^) 2 {l - 2^ 2 ) + 29^1^2 - Hi^d^xd^] (190) 

The coupling g a > flows to zero at large distances. On the other hand, we can absorb it by rescaling all fields so the 
action reads 

S=^[d 2 x [(^0)2(1 - 23^1^2) + 28^0^2 + Ag^^d^d^] (191) 

where now <f> has a different radius, <f> = 4> + "^f=- We see that as g a — » all interaction terms disappear and we get 
a free boson (f> together with a pair of free symplectic fermions ■01,2- Moreover the radius of compactification goes to 
infinity in that limit, so the boson <f> appears as non compact. 

This holds in the true large distance limit. At intermediate scales, we can use the RG equation for the coupling 

M 



dg a m — 2n — 2 2 1 2 



-9% (192) 



c?log/ 2-7T '"' n" 

Writing more generally 

-agl (193) 



d log I 

we see that g a approaches its vanishing large distance value as 

— = \ + a log(Z/Z ) s a \og(l/l Q ) (194) 

9a 9 a 

Here, I is a characteristic dimensionless scale ratio, roughly of the order of the ratio of the scale at which one is 
observing the physics to the lattice cut-off. On the cylinder, I can be identified with the width in lattice units, I = N. 

In the limit of large I, we can estimate more precisely the contribution to the spectrum coming from the boson <p. 
Recall that for a free bosonic theory where the action is normalized as S = J (d^X) 2 and the field compactified as 
X = X + 2nR, the spectrum of dimensions [JJ|] is 

x e 2 m 2 R 2 

A + A= i^ + ^ (195) 
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Matching the normalization gives R 2 — Air / ' g a in our case, and thus we expect the scaled gaps coming from the bosonic 
degrees of freedom to read, at large distances : 

A + A= 4^1ogG//o) +m27ralOS ^ o) (196) 

In the limit I — > oo the dimensions become degenerate and the spectrum can be considered as a continuum starting 
above A + A = 0. To emphasize the latter point, consider the contribution to the partition function coming from the 
<j> degrees of freedom. The system is defined on a torus of periods I and I', with I' '/I = r. Denote q = cxp(2«7rr). 



Z<p = J_^ ? (e/_R+mfl/2) 2 /2 -{e/ R-mR/2) 2 /2 



' ' e,m 



R 1 / irR 2 \mT-m n2 



a/2 Vim r nn ^ , V 2 Im r 

' ' m,m 

fa R 1 

R ~ * 00 \/2 Vlni r 7777 

where 77 = 9 1 ^ 24 II^Li(l ~ 9™) = q 1 ^ 2i P{q)- Observe now that one can write : 



(197) 



00 o s V 

■ ' (198) 



Vim r 7777 Jo 

which can be interpreted as an integral over a continuum of critical exponents A = A = s 2 . In the partition function 
(|197p . R plays the role of the density of levels, and is proportional to the (diverging) size of the target space. 

Going back to the specific case of the OSP(2\2)/OSP(l\2) model, we set a — l/ir. We get the radius R 2 — 
41og(Z// ), and the contribution of the free boson <p to the spectrum is : 

* * 1 (m 2 ) 2 log(Z/7 ) (199) 



41og(Z/J ) 

with e-2, 77i 2 arbitrary integers. Suppose now we consider the sigma model with a combination of periodic and antiperi- 
odic boundary conditions for the symplectic fermions. These get untwisted, and add to the discrete spectrum of a 
free boson at the free fermion point, which is given by equation (| 1 95[) with a radius R = 1. The quantum numbers 
associated to the fermions will be denoted e± , mi . Thus we expect : 

A + A = (ei ) 2 + t™f + + (m 2 ) 2 log(Z// ) (200) 

We can now compare with formula (|159[) . It is not entirely clear how constrained the quantum numbers might be in 
the lattice realization we are considering. But observe that finite scaled gaps occur for m + — m_ = m and converge 
to 

m 2 1 1 
x= — + -g{n + +n-Y + -K{ 1 ,N){n + -n-f (201) 

When 7 tends to tt/2, equations (|158[) and (|171[) give the coupling constants : 

.9 = 1, K(n/2,N)~-\ (202) 

loffiV 
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The numerical results are compatible with A! = 1. This suggests the identifications : 

mi = n+ + ri- 

62 = n + — ri- 
m,2 — 2^ m + — m -) = 

Note that the quantum numbers are (weakly) correlated : the integers mi, e2 are such that (mi + e2) is even. 
5.3 Interpretation of the staggered six- vertex model for 7 < 71/ 2 

In preparation for the subsequent discussion, we will denote the effective (square) radius of the non compact boson by 
(R2) 2 (equal to 41og(//Z ) for 7 = ir/2) and the radius of the compact boson by (R\) 2 (equal to 1 for 7 = ir/2). Away 
from 7 = 7r/2, it would be reasonable to try and interpret our results in terms of a deformation of the supersphere 
sigma model. Several scenarios are a priori possible. From the numerical results, the compact direction clearly has a 
modified radius which becomes : 

(R 1 f=g = ^l (203) 

7T 

As for the non compact direction, a first possibility would be to consider a constant 01(7) in the foregoing discussion. 
Although it provides reasonable results, much better fits are obtained with a radius going like log(//?o) for 7 7^ 7r /2; 
specifically, we find : 

(R2) 2 = ^ 1 [\og{l/lo)] 2 (204) 

57 

Note that this becomes ill-defined as 7 tends to 7r/2, where there is crossover to a behaviour linear in \og(l/lo). 

The most naive interpretation of the corresponding target space would be a torus with one period diverging like 
log(//Zo) in the scaling limit. This is reminiscent of results on the sausage model [20] which is a deformation of the 
usual sphere sigma model. There however, the theory instead of being massless in the IR is massless in the UV, while 
the target space in that limit is asymptotically a cigar. The dependence of the cigar dimensions on the RG scale and 
the anisotropy parameter are however reminiscent of ours; in particular the "long dimension" goes as the logarithm 
of the RG coordinate in both cases. 

Another, more suggestive interpretation, can be obtained if we consider the partition function of our model on 
the torus. Set 7 = ir/t so (R2) 2 = ^p[log(Z/Z )] 2 . Using the continuum representation (|19T|) and reabsorbing the 
t — 2 prefactor that comes from the dependence of R2 upon t into the continuously varying exponents we obtain the 
partition function : 

z J^m f fi,,^^ (205) 



e,m— — 00 



where the proportionnality constant is a presumably non universal quantity, independent of t, equal appproximately 

24 
47*^2) 



to 1/%/TO. Observe now that t^jp- — Trrzm = 2 so the conformal weights can be written as well, with respect to 



2(«+l) 
t-2 ' 



s 2 + j (m =F te) 2 

h = i + v ± ' (206) 

t-2 At K ' 
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This coincides with the contribution of the continuous representations to the spectrum of the Euclidian black hole 
CFT SL{2 1 R)/U{1) coset model [21]. 

It might be useful here to give some quick reminders for the spectrum of this model. The central charge for level 
k is : 

2(fc + l) , 9fm 
BH = k-2 ( ' 

Normalizable operators come in two kinds. There is the ones (delta functions normalizable) associated with the 
principal continuous series j ' = — A + is with conformal weight : 

, + , (n±kw) 2 s 2 + \ (n±kw) 2 

h --n—Y + 4 fc ~j—2 + 4 fc (208) 

and the ones coming from the discrete series. They have j € — 5] together with rules relating allowed values to 

w,n. The net result is that these fields all have dimensions larger than the bottom of (|208|) . Note that the identity 
field h = is not among the normalizable states, which is consistent with the fact that we do not observe (after the 
obvious identification k = t) cbh in the lattice model but cbh — 24 4 ^ 1 _ 2 - ) = 2. 

Of course the partition function of the Euclidian black hole theory should naively be infinite since it involves infinite 
dimensional representations of SL(2). The introduction of a Liouville wall at finite distance in the target space [22] 
gives a density of states p(s) oc loge to leading order, which agrees with our results provided we identify e = l/lo, and 
thus the size of the system (it would certainly be interesting to investigate the subleading behaviour, which depends 
on s, in the lattice data). 

Finally, we note that in [22] , the sum (4.17) implies some combinatorial contraints on the descendents, which we 
have not studied in the lattice model. 

Note that further twisting and reduction of the vertex model gives rise to parafermionic theories, which are 
themselves cosets SU(2)/U(1). It is not clear how this might be related to the identification of the untwisted vertex 
model with SL{2, R)/U(l). 



6 Geometrical interpretation of the critical exponents 

For the OSP(m\2n) models in their Goldstone phases, the continuous spectrum of critical exponents has its origin 
(like in the case of a pure non compact boson) in the existence of infinitely many operators with vanishing scaling 
dimension (the powers of <fi m the case of the non compact boson) , which in itself is a consequence of the spontaneously 
broken symmetry [7J. An obvious question is whether a similar interpretation exists for our model. 

Some insight can be gained by considering the limit 7 — + tt/2 which is related to the OSP{2\2) model. In the 
latter, exponents of the order parameter are related with geometrical correlations of the degrees of freedom carrying 
lines, and therefore it is tempting to ask whether this might hold away from the limit point as well. Before discussing 
this, it is important to notice that the staggered six- vertex model at 7 = n/2 (and the equivalent 38- vertex model) 
with periodic boundary conditions correspond, strictly speaking, to the model of [7J [TU] where the OSP symmetry is 
broken by the boundary. This is the reason why the spectrum of conformal weights contains a compact boson with 
(i?i) 2 = 1, and thus non trivial, finite exponents. Within the geometrical interpretation to be discussed below, this 
will correspond to the existence of some correlators having non trivial weights, while others do behave as in [7J. 

Let us now be more specific. As shown in section 12. 3[ summing on 2 x 2 vertex blocks and choosing the right 
basis, the staggered six- vertex model is equivalent to the 38- vertex model. Any lattice configuration within this model 
is completely covered by polygons of three different types : oriented lines, bare thin lines and bare thick lines. The 
particular point 7 = n/2, as well as our discussion on Goldstone phases, suggests to consider geometrical correlations 
associated with these lines. 
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6.1 Correlation functions with no thin or thick line 

Consider first the watermelon correlation function (as illustrated on the figure 131 [) consisting of a forced even number 
2r of positively-oriented lines (that is, we force 2r such lines to propagate through the system, on top of fluctuating 
numbers of positively and negatively oriented lies in equal number and bare lines). The critical exponent of this 
correlation function is given by the lowest-energy state in the sector with fixed total spin S — 2r and (7 = 1. 
According to equation Q159p . this state is v?(r,r).(o,o)- According to equation (|105[) . this is an eigenstate of the operator 
C, with eigenvalue one. 

6.2 Correlation functions with one thin or thick line 

Consider the correlation function consisting of a forced odd number 2r — 1 of positively-oriented lines, along with one 
bare (thin or thick) line. The critical exponent of the correlation function with a thin (resp. thick) line is given by the 
lowest-energy state in the sector with fixed total spin S = 2r — 1 and (7 = 1 (resp. C = —1). Both states tffc. ^ s 
have the physical exponent £( r ,r— i),(o,o)- Thus, the correlation functions with one thin or thick line have the same 
exponent £( r ,r-i),(o,o)- 

6.3 Correlation functions with one thin line and one thick line 

Consider the correlation function consisting of a forced even number 2r of positively-oriented lines, along with one 
thin line and one thick line. The critical exponent of this correlation function is given by the lowest-energy state in the 
sector with fixed total spin S — 2r and C = — 1. This state is <f7 r+1 r _n ( y with physical exponent £(r+i,r-i),(o.o)- 
In the particular case r = 0, this exponent becomes : 

X(i,-x U0 , )=K^,N) (209) 

This exponent vanishes in the thermodynamic limit. So the two- leg correlation function consisting of one thin line 
and one thick line belongs to the continuous subspectrum associated to the ground state. 
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6.4 Interpretation 

The physical interpretation we suggest for the continuum limit of the 38- vertex model is thus the following. The arrow 
degrees of freedom can be treated like the usual domain boundaries for a RSOS model which renormalizes, at large 
distances, to a compactified free bosonic field with radius R\. A correlation function involving a certain number S of 
positively oriented lines corresponds for this bosonic field to a magnetic or vortex operator, whose physical exponent 
is given by x — gS 2 /A. 

Meanwhile, correlators involving in addition thin and thick lines as well have, in the thermodynamic limit, exponents 
entirely determined from the contribution of the arrow degrees of freedom^ Thin and thick lines correspond thus to 
operators with vanishing exponents and (presumably) logarithmic correlators, and behave similarly to the crossing 
lines in the models of [TJ. 



7 Conclusion, open problems 

First and foremost, we believe this work puts on firm ground the existence of a continuous spectrum of critical 
exponents in a model with a finite number of lattice degrees of freedom per site (link), justifying fully the conclusions 

of mi]. 

It is truly remarkable that a proper staggering of the simple six-vertex model could give rise to such interesting 
behaviour: obvious directions for future work are plenty, and include an analytic derivation of the coupling constant K, 
a better understanding of the relationship with SL(2, R)/U(l), and of the effect of the various twist and truncations 
necessary to produce the Potts and RSOS versions of this model. 

It should also be possible to generalize the problem to some higher spin version. Some comments on Bethe equations 
are in order here. Recall that the usual source term for these equations reads, in the case of spin 1/2 : 



/ sinh | (a — i'y) \ N (210) 
\sinh | (a + i'y) J 



One might think of changing it through real heterogeneities ±A, leading to equations which have been used a great 
deal in the study of massive deformations 23J : 

/ • i 1 / . ■ \ \ N/2 , . . -i , . . „ . N/2 
( sinh i ( a — A — i'y) \ ( sinli ^ ( a + A — i'y )\ 



\ sinh i (a — A + i'y) J \ sinh i (a + A + ij) 

One could also think of changing it through imaginary heterogeneities. The simplest case would correspond to adding 
"string heterogeneities", i.e., for the simplest case of the two string, formally A = i'y. Clearly however half the modified 
source terms cancel out, leaving : 

/ sinh i(a — 217) \ 



Vsinhi(a + 2i 7 )J 



A subtle remark is in order here. Recall the conserved quantities of the 38-vertex model : the total spin and the value of operator C. 
In geometrical terms, this means that the total arrow flow is conserved, whereas only the parity of the number of thick lines is conserved. 
To define a correlation function with more than one thick lines within the 38-vertex partition function, it is necessary to define the transfer 
matrix on a non-local Hilbert space (so that the transfer matrix keeps track of the connectivity of the thick lines). In this framework, three 
interpretations are possible for the vertex . The behaviour of the correlation functions is likely to depend on the "splitting" of this 
vertex into the three possible connectivities, though we believe that the behaviour is universal, in agreement with this conclusion. We did 
not enter into these technicalities here, and simply described the correlation functions which are not affected by this problem. 
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which is nothing but the source term for a spin one chain. In general, heterogeneities of the type p-string will lead 
to the equations for the spin p/2 chain. The other natural possibility would be to add heterogeneities of the anti- 
string type, ie A = in. Up to a shift a — > a + in, this produces however the initial equations, and is well known to 
simply change the Hamiltonian from ferromagnetic to antiferromagnetic. The possibility we encountered in this paper 
consists in adding heterogeneitiesright in the middle of the "physical strip", at A = in/2. Note that the higher spin 
generalization is obvious by changing 17 to 2sij everywhere. 

It is also fascinating that the model should interpolate between OSP(2/2) and 5*0(4) symmetries when 7 goes 
from ir/2 to zero. This suggests the existence of a possible quantum group symmetry all along the line, which we have 
unfortunately not yet been able to identify. 
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